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DESCRIPTION 
A METHOD FOR SOLVING TRANSIENT SOLUTION AND 
DYNAMICS IN FILM BLOWING PROCESS 

5 TECHNICAL FIELD 

The present invention concerns a dynamic scheme for 
the film-blowing process and a method for solving transient 
solutions for the process. More specifically, the present 
invention solves the governing equations that take into 

10 consideration the viscoelasticity and cooling 

characteristics of film, and then through coordinate 
transformation, the invention transforms a free-end-point 
problem into a fixed-end-point problem. Then, by 

introducing Newton's method along with the OCFE (Orthogonal 

15 Collocation on Finite Elements) method, a new method for 
solving transient solutions is formulated. 

BACKGROUND ART 

The film-blowing process is a typical bi-axial 
20 extensional deformation process that produces oriented film 
by stretching and cooling polymer melts continuously 
extruded from an annular die in both axial and 
circumferential directions simultaneously, as shown in Fig. 
1. The axial extension is imposed by the drawing force of 
25 the nip rolls whereas the circumferential extension is 
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imposed by the air pressure inside the bubble. This film- 
blowing is similar to fiber spinning and film casting in 
engendering the extensional deformation of the material, 
yet salient in causing a biaxial extension. By manipulating 
5 the two important parameters of the system, i.e., the 
drawdown ratio (the ratio of the film velocities at the die 
exit and nip rolls) and the blowup ratio of the bubble 
between the die exit and the maximum bubble radius point, 
the process can be controlled as desired with respect to 

10 the process and the film product. 

Over the past four decades, many theoretical and 
experimental studies have been conducted on this important 
process. Among the major research results, the most 
comprehensive stability analysis was first carried out by 

15 Cain and Denn [Polym. Eng. Sci. 28: 1527, 1988] and then 
followed by Yoon and Park [Int. Polym. Proc. 14: 342, 1999]. 
Along with many interesting stability findings, the draw 
resonance instability, a self -sustained limit cycle-type 
supercritical Hopf bifurcation, has been well documented in 

20 these studies. 

While the basic understanding of the process in 
terms of steady state operations and linear stability has 
been greatly advanced by all these efforts, there still 
remains the need for a transient solution for the process 
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to reveal its nonlinear dynamics and nonlinear stability, 
which are acutely warranted for devising any systematic 
strategies for process stabilization and optimization. 
Unlike steady-state solutions that are relatively easy to 
5 obtain, the transient solutions of the governing equations 
of the process have long eluded theoretical pursuit, mainly 
due to the complex nonlinear nature of the partial 
differential equations and the boundary conditions. 

Due to the characteristics of a non-Newtonian fluid 

10 in terms of complex structure and difficulty in being 
interpreted, which is used as a subject for the 
rheologically-governed process, there is a variety in flow 
characteristics and instability. In particular, the 
representative instability, a draw resonance phenomenon 

15 that takes place during the extension deformation process, 
impedes the productivity of the process. Therefore, the 
process should be interpreted from a dynamic perspective, 
and control and design technology based on the nonlinear 
theory needs to be urgently developed more than anything 

20 else in order to overcome the instability and improve the 
productivity of the rheological process. Despite the 
significance of the film blowing process in industrial use, 
which produces films with wider width through the biaxial 
extension caused by the velocity difference and the 
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pressure difference, among all the extension deformation 
processes, the results of a transient solution for the 
process have not yet been reported due to the highly 
complicated nonlinearity of the governing equation, in 
5 comparison with the other processes. So far, neither 
nonlinear stability analysis nor a transient solution for 
the non-isothermal governing equations, which explain the 
cooling of the film-blowing process, have been reported 
anywhere in the world. 

10 The inventors employed the Phan Thien-Tanner (PTT) 

constitutive model, known for its capability to accurately 
portray the extensional flows of viscoelastic polymers, and 
an energy equation where the cooling characteristics of the 
film is taken into account for the system, thereby enabling 

15 them to formulate a transient solution for the non- 
isothermal film-blowing process and analyze the stability 
of the nonlinear system. 

Disclosure of the Invention 
20 TECHNICAL PROBLEMS 

The purpose of the present invention is to solve the 
governing equations in the form of partial differential 
equations for the film-blowing process with its 
viscoelasticity and cooling characteristics in mind. Then, 
25 through coordinate transformation, we transformed the free- 
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end-point problem into a fixed-end-point problem, and 
finally, by using numerical schemes such as Newton' s method 
and OCFE (Orthogonal Collocation on Finite Elements), we 
provide a new method for yielding transient solutions. 

5 

TECHNICAL SOLUTIONS 

The abovementioned objective of the invention is to 
yield a governing equation for the film-blowing process, in 
which viscoelasticity and cooling are considered, and 

10 effectively simulate the real process by effectively 
formulating a transient solution through methods such as 
Newton's method or the OCFE. It is expected that the 
results drawn from the calculation can be applied in 
developing a device for the optimization and stabilization 

15 of the film blowing process, establishing optimal operating 
conditions, and developing a polymer material. 

The present invention is characterized by four 
stages. The first stage is the solving of governing 
equations. The next stage is the yielding of transient 

20 solutions through the use of coordinate transformation and 
numerical methods such as Newton's method and OCFE. The 
third stage is the comparison stage wherein the results 
from the calculations are compared with the actual results 
from the experiments of the process. In the final stage, 
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strategies for the optimization and stabilization of the 
process are drawn up. 



The present invention solves the governing equations 



transformed along the temperature-time axis and then, by 
using Newton's method and the OCFE (Orthogonal Collocation 
on Finite Element) method, we provide the method for 
solving transient solutions for the non-isothermal film- 
blowing process: 

Equation : 



in the manner shown below. 



Then, the coordinates are 
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Here, 




Axial direction: 




• • • (2) 
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Circumferential direction : 



. . . (3) 



Constitutive Equation : 



. (4) 



10 Here, 

jr=exp[cZfetrr], £= V»~£a 2ZM Vv+ ViA 2*^= -^2-. 



Energy equation: 

15 . . . ( 5 ) 

Here, 
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Boundary conditions : 
H« ijjw |*s £|ss l^T 8889 !^ 



at z=0 . . . (6a) 



J?IL =4- >,.^..v.v^ ,....».,.„ , . .J . i l?rES ^ BBS II , „■„ ■■ i mm jp«- 0wm ff- 



at z=z F . . . ( 6b) 



In the above equations, r denotes the dimensionless 
bubble radius, w the dimensionless film thickness, v the 
dimensionless fluid velocity, t the dimensionless time, z 

10 the dimensionless distance coordinate, AP the air pressure 
difference between inside and outside the bubble, B the 
dimensionless pressure drop, A the air amount inside the 
bubble, P a the atmospheric pressure, T z the dimensionless 
axial tension, C gr the gravity coefficient, T drag the 

15 aerodynamic drag, a SU rf the surface tension, 9 the 
dimensionless film temperature, i the dimensionless stress 
tensor, D the dimensionless train rate tensor, £ and £ the 
PTT model parameters, De the Deborah number, 9o the zero- 
shear viscosity, K the dimensionless activation energy, U 

20 the dimensionless heat transfer coefficient, E the 
dimensionless radiation coefficient, k a ± r the thermal 
conductivity of cooling air, p a ± r the density of cooling air, 
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rj a ± r the viscosity of cooling air, v c dimensionless cooling 
air velocity , a and (5 parameters of heat transfer 
coefficient relation, 9 C the dimensionless cooling-air 
temperature, 0«, the dimensionless ambient temperature, e m 
5 the emissivity, a SB the Stefan-Bolt zmamn constant, p the 
density, C p the heat capacity, and D R the drawdown ratio. 

However, the assumption was made that no deformation 
occurred in the film past the freezeline at the boundary 
conditions. Overbars denote the dimensional variables. 

10 Subscripts 0, F and L denote the die exit, the freezeline 
conditions and the nip roll conditions, respectively. 
Subscripts 1, 2 and 3 denote the flow direction, normal 
direction and circumferential direction, respectively. 

In the present invention, while in the process of 

15 yielding numerical solutions for the isothermal film- 
blowing process, we employed a coordinate transformation to 
make time and temperature as new independent variables in 
lieu of the original time and distance. Through this 
transformation, the moving freezeline height can be handled 

20 effectively because the boundary conditions are clearly set. 
Also, in the numerical analysis methods, the Newton's 
method and OCFE were employed. 

In the non-isothermal process model used for the 
present invention, a total of three multiplicities were 
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discovered, and this matches the experimental results. In 
addition, after analysis of the three points of 
stabilization where multiplicities exist, we discovered 
that at points in the bubble with the smallest and the 
5 largest radius, the disturbances introduced into the system 
disappear after some time. However, at points in the 
middle, the disturbances increase and draw resonance occurs. 
This also matches the experimental results. Especially in 
the case of draw resonance, the transient solutions for the 
10 amplitude and period of the bubble radius can be exactly 
predicted. 

Hereafter, practical examples will be used to 
explain in more detail the specific methods of the present 
invention. However, the applications of the present 
15 invention are not limited to only these examples. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows the film-blowing process. 

Fig. 2 is a graph that shows the optimum number for 
20 each factor along with the optimum number for the 
collocation point by making use of the transient response 
indications. Inaccurate results: NE=4, NP=5 (-•-); NE=5, 
NP=4 (-••-); Accurate results: NE=5, NP=5(-); NE=5, NP=6 ( — 
); NE=6, NP=5 ( • ■ •) . (a) shows the case during small time, 
25 (b) during large time, and (c) when NE=5 and when NP=5. 
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Fig. 3 represents the typical variations of the 
bubble radius and the freezeline height along the flow 
direction during one period of the draw resonance 
oscillation, when the sustained periodicity of the draw 
5 resonance is fully developed. 

Under the same condition of Fig. 3, Fig. 4 shows the 
dimenionless bubble radius changes during one period of 
oscillations plotted (a) against the dimensionless distance 
from the die exitto the freezeline, z, and (b) against the 
10 transformed temperature coordinate (zeta). 

Fig. 5 shows the multiplicities of the non- 
isothermal process obtained by observing the intersections 
of the straight line of a constant drawdown ratio (D R =35) 
and the curve of a constant air pressure (B=0.37). (a) 
15 shows the simulations results, (b) shows the stability at 
point H, (c) shows the numerical scheme results at point H 
along with the draw resonance that appeared during the 
experiment, and finally (d) shows the stabilization at 
point L. 

20 

BEST MODE 

Example 1 : Solution for the governing equations of 
the film-blowing process and yielding of transient 
solutions by using numerical methods 
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Trying to solve the new governing equations while 



considering the non-isothermal characteristics of the film- 
blowing process and then yielding transient solutions by 
using several numerical analysis methods is vital in the 
5 theoretical study concerning the stabilization of the 
process . 



isothermal film-blowing of PTT fluids, based on the seminal 
work of Pearson and Petrie {J. Fluid Mech. 40(1970) 1 and 
10 42(1970) 609), who established the first modeling equations 
and the standard for all ensuring research efforts, are as 
follows : 

Equations : 



The dimensionless governing equations of the non- 




15 



. . . (1) 



Here, 




Axial direction: 




20 



. . . (2) 



Here, 
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Circumferential direction : 



— — — — , 

... (3) 



Constitutive Equation: 

j 



(4) 



10 Here, 



Energy Equation: 

j$£ i iy3 rjr J? 1 

Vi+Xar/azf ** ™ * * 

15 . . . ( 5 ) 

Here, 
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5 



Boundary Conditions : 

^ p ^S ^ S£ j^^iS^f^ 

at z = 0 ... (6a) 

— T + — ^=s=^==r w» 0. jBssssssssr *■* 23W $U 

at »z vT+{i^F v^Ri^F * ' 

at z=z F . . . ( 6b) 



In the above equations, r denotes the dimensionless 
bubble radius, w the dimensionless film thickness, v the 

10 dimensionless fluid velocity, t the dimensionless time, z 
the dimensionless distance coordinate, AP the air pressure 
difference between inside and outside the bubble, B the 
dimensionless pressure drop, A the air amount inside the 
bubble, P a the atmospheric pressure, T z the dimensionless 

15 axial tension, C gr the gravity coefficient, T^rag the 
aerodynamic drag, a surf the surface tension, 0 the 
dimensionless film temperature, t the dimensionless stress 
tensor, D the dimensionless train rate tensor, e and £ the 
PTT model parameters, De the Deborah number, 0 O the zero- 

20 shear viscosity, K the dimensionless activation energy, U 
the dimensionless heat transfer coefficient, E the 
dimensionless radiation coefficient, k air the thermal 
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conductivity of cooling air, p a ir the density of cooling air, 
rj a ±r the viscosity of cooling air, v c dimensionless cooling 
air velocity , a and (3 parameters of heat transfer 
coefficient relation, 9 C the dimensionless cooling-air 
5 temperature, 9~ the dimensionless ambient temperature, e m 
the emissivity, o S b the Stefan-Bolt zmamn constant, p the 
density, C p the heat capacity, and D R the drawdown ratio. 

However, the assumption was made that no deformation 
occurred in the film past the freezeline at the boundary 

10 conditions. Overbars denote the dimensional variables. 
Subscripts 0, F and L denote the die exit, the freezeline 
conditions and the nip roll conditions, respectively. 
Subscripts 1, 2 and 3 denote the flow direction, normal 
direction and circumferential direction, respectively. 

15 Several assumptions have been incorporated in the 

above model: 

First, the thin film approximation that all state 
variables depend on the time and z-coordinate, simplifies 
the system into a one-dimensional model. 
20 Second, the bubble is axisymmetric, excluding 

possible helical instability. 

Third, the secondary forces acting on the film, such 
as inertia, gravity, air-drag and surface tension are 
neglected. 
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Fourth, the crystallization kinetics of polymer 
melts are not included here. 

Finally, the origin of the z-coordinate is chosen at 
the point of extrudate swell, meaning the deformation of 
5 polymer melts inside the die being lumped into the initial 
conditions at z=0 . 

It has proven to be impossible to yield transient 
solutions in the above governing equations for the non- 
isothermal film-blowing process using conventional 
10 numerical schemes. Especially in cases where draw 

resonance instability exists in the process, a new and 
effective numerical scheme has to be devised in order to 
yield transient solutions. 

First, we tried a finite difference method (FDM) of 
15 successive iterations that involves solving each equation 
for one variable while the other state variables are 
assumed as known. Although this method has been successful 
in fiber spinning and film casting, it failed for the 
present invention mainly because of the existence of a 

20 nonlinear term in the equations (i.e. 1 H~ (S^/ote) ) which 
stems from the fact that the fluid velocity is in the film 
direction, not in the machine direction. 

Next, we applied the Newton's method with FDM to 
simultaneously solve the equations for all dependent 
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variables. However, this method entails an extremely long 
computation time in the order of weeks, if computations are 
ever possible, due to the full matrix calculations, thus 
rendering itself unworkable for all practical purposes, 
5 Finally, we introduced an orthogonal collocation 

method on the finite elements of z-coordinate (OCFE) . 
Employing a minimum number of finite elements (NE) and a 
minimum number of collocation points (NP) within each 
element to guarantee accurate transient solutions with 

10 manageable computation times, both of which turned out to 
be five in the present invention (Fig. 2a and 2b), we 
finally succeeded in devising a numerical scheme for 
generating transient solutions for the process even during 
the instability of draw resonance. Analytically-derived 

15 expressions for each element in the Jacobian matrix further 
facilitate the solution procedure with much ease. For the 
transient simulation, an implicit second-order backward 
scheme in time-derivative terms was used to enhance 
numerical robustness. Fig. 2c shows a typical example of 

20 the time convergence of a transient solution in draw 
resonance . 

With the OCFE, we also introduced several important 
modeling ideas for a more accurate description of the 
system. First, to handle the moving freezeline height, a 
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coordinate transformation was employed to make time and 
temperature as new independent variables in lieu of the 
original time and distance. This transformation essentially 
converted the free-end-point problem into a computationally 
5 amenable fixed-end-point one. The following coordinate 
transformation was applied: 

<= — — 

Here, the new independent variable £ becomes 0 at 
the die exit and 1 at freezeline height. By applying the 
10 above transformation to the governing equations, a new 
(t, O coordinate replaces the (t,z) coordinate: 

(MA = (M\ /j^r 1 

\3z$ t \6CU\OCh 

... (8) 

\ at } m \ at / c \ &c | t \ &cl t \&t) € 

... (9) 

Here, f represents all state variables. 

15 Second, instead of the so-called cylindrical 

approximation in calculating the amount of air pressure 
inside the bubble as used by others, in the present 
invention, the actual shape of the bubble is traced when 
calculating the real bubble volume. This allows the exact 

20 temporal shape of the propagating bubble disturbances to be 
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captured in the simulation during the oscillating 
instability . 

Fig. 3 shows the comparison of simulation data in 
draw resonance with a real experimental case. To our 
5 knowledge, this demonstration of transient behavior is the 
first in the literature. In view of the assumptions 
incorporated in the modeling of the highly nonlinear 
dynamical process of film-blowing, the closeness of the 
simulation results to real observations is considered as a 

10 modeling and numerical breakthrough. To clearly depict the 
transient behavior of the state variables in this draw 
resonance, the dimensionless bubble radius during one 
period of the oscillation is plotted in Fig. 4a against the 
dimensionless distance from the die exit to the freezeline 

15 height, i.e., the original independent variable (z) , and 
also plotted in Fig. 4b against the transformed 
dimensionless temperature coordinate (Or i.e., the new 
independent variable, which always has the same unity value 
at the freezeline point. 

20 Fig. 5 exhibits an interesting case where three 

experimentally observed steady states were simulated quite 
closely, attesting to the usefulness and robustness of the 
simulation model. The three steady states in these 
particular cases were determined in the stability diagram 
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by the intersections of the straight line of a constant 
drawdown ratio (which has a fixed slope of 1/35 because the 
drawdown ratio ( DR = 35 here) is, by definition, equal to 
the ratio of the thickness reduction (TR) and the blowup 
5 ratio of the bubble (BUR) ) and the curves of a constant air 
pressure inside the bubble (B = 0.37 here). The stability 
diagram in Fig. 5a was obtained using linear stability 
analysis. Among these three steady states, only the middle 
one (Fig. 5c) turns out to be unstable, exhibiting draw 

10 resonance, whereas the other two steady states, the upper 
and lower BUR steady states (Fig. 5b and d) are stable. In 
non-isothermal film-blowing, not only the bubble radius, 
but also the other state variables such as film thickness, 
bubble air pressure and freezeline height, all oscillate 

15 with time during draw resonance instability. The typical 
oscillation results of the bubble radius at the freezeline 
are shown in Fig. 5c, exhibiting an excellent agreement 
between the off-line film experimental data and the 
theoretical on-line simulation data. 

20 The utility of these transient solutions for the 

film-blowing process is rather far-reaching in both the 
analysis and synthesis of the system. First, it enables us 
to confirm the same draw resonance criterion previously 
developed in fiber spinning and film casting based on the 
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traveling times of kinematic waves, also applying to film- 
blowing. Second, the sensitivity analysis in assessing the 
effects of process conditions such as cooling, 
viscoelasticity of input polymers, and the air 
5 amount/pressure inside the bubble, on the behavior of the 
system can be easily performed with transient solutions as 
in other extensional deformation processes. Third, taking 
advantages of the two utilities mentioned above, we will be 
able to develop strategies for finding the optimal 

10 conditions for cooling, polymer viscoelasticity, air 
pressure/amount, freezeline height, etc., leading to 
enhanced productivity and film quality. Fourth, the 
transient solutions can also be applied in developing the 
apparatuses necessary for the optimization and 

15 stabilization of the process and can also prove useful in 
developing polymer materials. 



ADVANTAGEOUS EFFECTS 

As we have investigated through the practical 
20 example shown above, the present invention concerns a 
method for the dynamic scheming and yielding of transient 
solutions in the film-blowing process. By using the 
Newton's method and OCFE (Orthogonal Collocation on Finite 
Element), we were able to achieve transient solutions for 
25 the non-isothermal film blowing process. These solutions 
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proved to be of great use when analyzing the nonlinear 
stability and nonlinear dynamics of the process. Also, 
through experiments, the numerical solutions that were 
obtained theoretically were verified to be useful. By 
5 taking into consideration the dynamics of the process that 
show nonlinear motion, the transient solutions may be 
applied in the optimum design of the process and in 
nonlinear control. We can also develop the apparatuses 
necessary for the optimization and stabilization of the 
10 film-blowing process, thereby realizing high-productivity 
and high-quality products. Therefore, the present 

invention is highly valuable in processes such as film, 
coating and flat display wherein the transformation of the 
film is important. 
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